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Preface 


A  three-component  magnetometer  study  was  undertaken  as  part  of  the  Work 
Unit  “Innovative  Geophysical  Technologies  for  Enhanced  Buried  UXO 
Discrimination”  (AF25, 6.2),  funded  by  the  U.S.  Army  Engineer  Research  and 
Development  Center  (ERDC),  under  die  Environmental  Quality  Technology  (EQT) 
Program.  The  data  were  acquired  during  the  period  January  through  April  2000  by 
personnel  of  the  Geotechnical  Laboratory  (GL)  and  Information  Technology 
Laboratory  (ITL),  ERDC,  Vicksburg,  MS.  This  was  a  preliminary  investigation  to 
determine  if  the  additional  information  provided  by  three-component  magnetometer 
sensors  can  improve  the  discrimination  and  identification  capabilities  currently 
achieved  using  total  magnetic  field  magnetometers.  Dr.  M.  John  Cullinane  was 
Program  Manager,  USAERDC,  for  the  EQT  Program. 

This  report  was  prepared  by  Dr.  Janet  E.  Simms,  Geotechnical,  Earthquake,  and 
Geosciences  Division  (GEGD),  GL.  Field  data  were  collected  by  Dr.  Simms  and 
Mr.  Lewis  B.  Smithhart,  Operations  Branch,  ITL.  Geophysical  data  analysis  was 
performed  by  Dr.  Simms.  The  work  was  performed  under  the  direct  supervision  of 
Mr.  Donald  Yule,  Acting  Chief,  Geotechnical,  Earthquake,  and  Geophysics  Branch, 
GEGD,  and  the  general  supervision  of  Drs.  M.  E.  Hynes,  Acting  Chiefs  GEGD,  and 
Michael  J.  O’Connor,  Director,  GL. 

At  the  time  of  publication  of  this  report.  Director  of  ERDC  was  Dr.  James  R. 
Houston.  Commander  was  COL  James  S.  Weller,  EN. 


The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication , 
or  promotional  purposes .  Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of the  use  of  such  commercial  products. 


Conversion  Factors, 
Non-SI  to  SI  Units  of 
Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units  as 
follows: 


1  Introduction 


The  objective  of  the  Work  unit  “Innovative  Geophysical  Technologies  for 
Fnhamwt  Buried  UXO  Discrimination”  is  to  develop  advanced  geophysical 
technologies  to  enhance  the  ability  to  discriminate  buried  UXO  in  a  wide  range  of 
environmental  and  geophysical  conditions.  The  geophysical  methods  under  study 
include  multi-component  magnetic,  gravimetric,  and  time  domain  electromagnetics. 
This  report  details  the  initial  findings  using  multi-component  magnetic  sensors. 

With  the  closing  of  numerous  installations  and  increased  environmental 
awareness,  a  primary  thrust  of  the  Army  is  the  cleanup  of  unexploded  ordnance 
(UXO)  from  military-owned  land  and  transfer  of  this  land  to  the  civilian  sector. 
Geophysical  methods  are  commonly  employed  for  locating  UXO.  Of  these 
magnetometers  (total  field  and  gradient)  and  “simple”  time  domain 
electromagnetic  (TDEM)  instruments  are  routinely  used.  Both  types  of  instruments 
are  excellent  for  the  detection  of  potential  UXO,  which  are  usually  comprised  of 
both  ferrous  and  non-ferrous  metals.  The  ubiquity  of  the  UXO  metallic  make-up 
and  the  nature  of  the  geophysical  instrument  measurement  makes  it  very  difficult  to 
^icrriminate  potentially  dangerous  UXO  from  non-threatening  exploded  ordnance 
fragments  and  waste,  cultural  objects,  and  anomalous  geologic  sources.  A  large 
percentage  of  the  cost  (approximately  75%)  of  UXO  cleanup  efforts  is  attributed  to 
the  detection  and  excavation  of  “false  alarm”  anomalies,  i.e.,  anomalies  caused  by 
buried  ordnance  debris,  other  metallic  objects,  and  naturally  occurring  soil  features 
(USAEC  19  JPG  reference).  Therefore  it  does  not  suffice  to  just  detect  an  anomaly, 
it  is  necessary  to  discriminate  and,  ultimately,  identify  UXO. 

The  term  “simple”  used  above  to  describe  the  TDEM  instruments  refers  to 
those  instruments  that  measure  only  one  or  two  time  gates  within  the  transient  decay 
window.  The  consensus  among  many  UXO  practitioners  is  that  the  total  magnetic 
field,  magnetic  gradient,  and  “simple”  TDEM  measurements  do  not  provide 
ciiffiriwit  information  to  adequately  discriminate  UXO  from  other  anomalies. 
Multi-component  magnetic  data,  multi-channel  TDEM  data,  or  multi-sensor  data 
likely  are  needed  to  obtain  a  more  accurate  description  of  the  anomaly  orientation 
(azimuth  and  dip),  size  and  depth  to  allow  discrimination  of  UXO  from  non-UXO 
anomalies  This  report  concentrates  on  multi-component  magnetic  data,  specifically 
three  component  magnetic  data  acquired  using  a  Barrington  Mag-03  sensor. 

The  majority  of  UXO  detection  work  that  utilizes  magnetometers  measures  the 
total  magnetic  field  or  gradient.  Those  surveys  that  do  use  a  multi-component 
magnetometer,  such  as  the  Schiebel  DIMADS™,  generally  do  not  evaluate  the 
individual  component  data  but  rather  calculate  the  total  magnetic  field  from  the 
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components  and  use  it  Much  research  has  been  done  in  evaluating  the  total 
magnetic  field  signature  ova:  UXO  (e.  g.  Altshuler  1996;  Barrow  et  al.  1996; 
Barrow  et  al.  1997;  Barrow  and  Nelson  2000;  Hollyer  et  al.  1997).  Early  work 
stressed  the  need  for  data  coverage  on  the  order  of  die  size  of  die  smallest  ordnance 
item  expected  to  be  detected.  These  results  were  the  impetus  for  designing  multi- 
array  sensor  configurations  (e.g.  STOWS,  MTADS).  Even  with  the  improvement  in 
data  coverage,  the  discrimination  of  UXO  from  ordnance  waste  and  other  man-made 
and  geologic  sources  has  lagged  detection  capabilities.  The  series  of  UXO 
Technology  Demonstrations  at  Jefferson  Proving  Ground  (USAEC  1994, 1995, 
1996, 1997;  Altshuler  et  al.  1995)  emphasize  this  issue.  Recent  advances  in  UXO 
discrimination  primarily  have  been  achieved  through  computer  analysis  of  total 
magnetic  field  data  alone  (modeling:  Altshuler  1996;  Barrow  2000;  Barrow  and 
Nelson  2000;  Barrow  et  al.  1997;  Butier  et  al.  1998;  McFee  and  Das  1990;  neural 
network  analysis:  Grimm  2000;  Hart  2000;  Millhouse  2000;  image  processing: 
Millhouse  2000)  or  through  joint  analysis  with  electromagnetic  induction  data 
(Collins  2000;  Grimm  2000;  Hart  2000;  Lavely  and  Grimm  1997).  The  present 
work  takes  a  different  approach  through  the  study  of  three-component  magnetic 
data.  Model  data  are  generated  and  compared  to  field  data  collected  under 
controlled  field  conditions.  The  data  obtained  for  the  three  components  are  analyzed 
to  determine  if  the  additional  information  provided  by  the  components  aids  in  the 
discrimination  of  buried  ordnance. 
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2  Theoretical  Magnetic  Models 


The  forward  magnetic  modeling  program  MAGMOD  was  used  to  generate  the 
response  of  a  105 -mm  projectile.  The  features  and  capabilities  of  MAGMOD  are 
described  in  Butler  et  al.  (1998)  and  Butler  et  al.  (2000).  MAGMOD  is  based  on 
the  theoretical  development  by  McFee  and  Das  (1990)  and  uses  (1)  a  prolate 
spheroid  and  (2)  a  multipole  expansion  solution  for  the  induced  magnetic  field 
external  to  the  spheroid  to  approximate  the  shape  of  a  projectile.  The  multipole 
expansion  includes  the  dipole  and  octupole  terms,  having  no  monopole  term  and  the 
quadrapole  term  is  zero  due  to  symmetry  of  the  spheroid.  Since  the  octupole  term 
falls  off  at  a  rate  of  1/r s,  the  effect  of  higher  order  terms  is  negligible.  Input  to  the 
program  includes:  the  earth’s  field  strength,  inclination  and  declination;  spheroid 
length,  diameter,  azimuth  and  dip;  x,  y  coordinates  and  depth  to  the  spheroid  center; 
relative  magnetic  permeability  of  model;  sensor  height  and  grid  dimensions. 

Figure  1  shows  the  model  total  magnetic  field  anomaly  generated  for  a  105 -mm 
projectile  as  the  azimuth  is  varied  from  0°  to  90°  (constant  dip)  (a)  and  when  the  dip 
is  varied  from  0°  to  90°  (constant  azimuth)  (b).  A  geomagnetic  field  having  a 
declination  of  0°  with  an  inclination  of  65°  was  used  in  constructing  the  models. 
Distance  between  the  sensor  and  center  of  the  spheroid  is  50cm.  The  data  were 
generated  over  a  5m  x  5m  grid  using  a  12.5cm  spacing.  The  image  plots  use  the 
game:  chatting  scale.  The  profile  plots  below  each  image  plot  represent  the  north- 
south  profile  across  the  center  of  the  image  plot  The  x  and  y  axes  of  the  profile 
plots  are  to  the  same  scale.  As  both  the  azimuth  and  dip  increase  toward  90°,  the 
anomaly  structure  changes  and  approaches  that  of  a  monopole.  Also  note  in  Figure 
la  that  the  magnetization  vector  (direction  of  dipole)  lags  that  of  the  azimuth 
(Altshuler  1996).  For  a  length  to  diameter  ratio  of  4,  which  is  the  aspect  ratio  of  a 
105-mm  projectile,  the  magnetization  vector  begins  to  lag  the  azimuth  at  azimuths 
greater  than  about  5°.  A  maximum  lag  of  22°  occurs  when  the  angle  between  the 
long  axis  of  the  projectile  and  the  geomagnetic  field  is  68°. 

A  uniformly  magnetized  sphere  (simple  induced  dipole)  is  the  simplest  model 
for  approximating  the  anomalous  response  of  a  UXO.  Figure  2  depicts  plots  of  a 
spherical  model  fit  to  the  profiles  in  Figure  1.  A  sphere  having  a  diameter  of  15cm 
was  arbitrarily  chosen  and  does  not  represent  an  equivalent  volume  to  the  105-mm 
projectile.  The  only  parameter  to  vary  using  a  spherical  model  is  depth  (z),  yet 
fairly  good  fits  to  the  spheroid  models  are  obtained,  especially  when  the  azimuth 
(with  dip  =  0°)  (a)  and  dip  (with  azimuth  =  0°)  (b)  are  90°.  However,  the  spherical 
model  does  not  provide  any  information  on  the  azimuth,  dip  or  shape  of  the  buried 
object  More  sophisticated  models,  such  as  a  full-field  solution  (Stratton  1941; 
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Altshuler  1996),  or  a  dipole  or  multipole  expansion  solution  for  a  spheroid,  or 
multi-component  data  are  needed  to  achieve  discrimination  and  identification  of 
UXO. 


The  total  magnetic  field  Ht  can  be  described  by  its  vector  components  Hx , 

Hy ,  Hz  and  its  orientation  with  respect  to  the  horizontal  and  vertical  planes.  By 

convention,  Hx  points  in  the  direction  of  true  north  (TN).  Referring  to  the  diagram 
below,  the  declination  D  is  the  angle  measured  in  the  horizontal  plane  between 

Hx  and  magnetic  north  (MN),  and  the  inclination  is  the  angle  the  total  magnetic 
field  vector  makes  with  the  horizontal  plane. 


X  TN 


Hx  =  Hh  cos  D 
Hy  =  Hh  sin  D 
Hz  =  Ht  sin  I 
Hh  =  Ht  cos  I 


Magnetic  component  space  plots  for  the  105 -mm  spheroid  discussed  above  are 
plotted  in  Figure  3.  The  component  data  correspond  to  the  north-south  profile 
acquired  across  the  center  of  the  image  plots  in  Figure  1,  where  each  point  on  the 
plots  in  Figure  3  correspond  to  the  (Hx,  Hz)  or  (Hx,  Hy)  values  for  a  point  on  the 
profile  lines.  Figure  3a  is  for  the  case  dip  =  0°  and  azimuth  varies  0°  to  90°.  The 
left  plot  represents  the  vertical  component  Hz  plotted  against  the  horizontal 
component  Hx,  whereas  the  right  plot  is  the  horizontal  component  Hy  versus  Hx. 
The  directional  components  are  oriented  based  on  a  right-hand  system  with  z 
positive  downward  (x  north,  y  east).  The  heart-shaped  Hz-Hx  curves  are  not 
symmetrical.  At  an  azimuth  of  0°  the  curve  lays  horizontal  (“heart”  on  its  side)  and 
as  the  azimuth  of  the  projectile  increases,  the  curve  rotates  toward  a  vertical 
position.  The  curves  in  die  horizontal  component  plot  Hy-Hx  are  quite  different 
depending  on  azimuth.  Due  in  part  to  symmetry  of  the  spheroid  and  the  fact  that  the 
geomagnetic  field  is  directed  in  the  x  (north)  direction,  no  variation  in  the 
y-component  is  observed  at  an  azimuth  of  0°  and  the  curve  is  a  horizontal  line.  In 
addition  to  these  two  factors,  when  the  spheroid  is  rotated  from  the  direction  of  the 
geomagnetic  field,  the  magnetization  vector  follows  but  lags  the  long  axis  of  the 
spheroid  up  to  a  maximum  angle,  at  which  point  the  direction  of  magnetization 
begins  to  move  back  toward  the  direction  of  the  geomagnetic  field.  At  an  azimuth 
of  90°  the  magnetization  vector  returns  parallel  to  the  geomagnetic  field  (Altshuler 
1996);  the  horizontal  line  for  a  90°  azimuth  in  the  Hy-Hx  plot  reflects  this 
phenomenon.  As  the  azimuth  of  the  spheroid  increases  from  0°,  a  magnetic 
influence  is  now  present  in  both  the  x-  and  y-components  and  a  “tornado-shaped” 
curve  forms  that  rotates  toward  vertical.  When  the  azimuth  is  held  constant  and  the 
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dip  of  the  projectile  is  varied,  Hz-Hx  plots  similar  to  those  when  the  azimuth  was 
varied  are  obtained  (Figure  3b).  The  heart  curve  is  relatively  horizontal  at  a  dip  of 
0°  and  rotates  to  vertical  as  the  dip  of  the  projectile  increases  to  90°.  The  Hy-Hx 
curves  are  a  horizontal  line  for  all  dips,  indicating  a  zero  y-component  which  is 
expected  for  an  azimuth  of  0°  and  a  north  oriented  geomagnetic  field. 
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3  Field  Acquisition  and  Data 
Review 


This  section  describes  how  the  field  data  were  obtained  and  presents  the  results 
of  the  three  magnetic  components  measured,  both  as  profile  plots  and  image  plots. 

Data  Acquisition 

Three  component  magnetic  data  were  acquired  under  controlled  field  conditions. 
A  Barrington  Mag-03MS  low  noise  sensor  having  a  measuring  range  of  ±100 
microtesla  was  used.  The  sensor  has  a  noise  error  of  6  picotesla  rms/VHz  and  an 
orthogonality  error  of  ±  0. 1%.  The  Mag-03  is  a  fluxgate  sensor  that  measures  three 
components  of  the  magnetic  field,  two  horizontal  and  one  vertical.  The  components 
are  oriented  in  a  right  hand  system  with  the  vertical  positive  downward.  The  survey 
grid  measured  4  meters  east-west  and  5  meters  north-south.  Data  were  collected  at 
a  height  of  23cm  above  the  ground  surface  with  a  horizontal  spatial  increment  of 
25cm.  Figure  4a  shows  the  measurement  stand  and  4b  the  sensor  holder.  A  bubble 
level  was  mounted  on  the  sensor  holder  and  a  wood  wedge  was  used  to  level  the 
sensor.  No  metal  was  used  in  the  construction  of  the  measurement  stand  nor  the 
sensor  holder.  Data  were  collected  over  a  105 -mm  projectile  (length  42cm,  diameter 
10.5cm)  in  a  horizontal  position  (dip  -  0°)  and  for  north-south  (nose  pointing  north) 
and  east-west  (nose  pointing  east)  orientations.  The  projectile  was  located  in  the 
center  of  the  grid  and  buried  at  a  depth  of  5.25cm,  measured  to  the  center  of  the 
projectile;  the  top  surface  of  the  105-mm  was  at  ground  surface.  A  base  station  also 
utilizing  a  Mag-03MS  sensor  was  used  to  correct  for  diurnal  variations.  The  data 
acquisition  module  for  the  Mag-03  sensor  allows  input  from  two  channels,  therefore 
a  base  station  reading  was  acquired  each  time  a  measurement  on  the  survey  grid  was 
collected. 

Data  Presentation 


The  component  data  plots  Hz-Hx  and  Hy-Hx  for  both  azimuthal  orientations 
are  presented  in  Figure  5.  The  data  are  from  the  north-south  profile  acquired  along 
the  centerline  of  the  grid.  The  model  curves  generated  using  MAGMOD  for  a  105- 
mm  projectile  and  using  the  local  earth  magnetic  field  parameters  (except  for  the 
declination  which  was  fixed  at  zero)  are  also  plotted  for  comparison.  Figure  5a 
shows  the  component  plots  for  an  azimuth  of  0°.  The  Hz-Hx  curves  are  similar  in 
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shape  for  both  the  model  and  field  data,  although  the  field  data  has  a  greater  range 
in  magnihiHft  The  model  curve  exhibits  an  inner  loop  that  is  caused  by  the 
influence  of  the  octupole  tom.  The  octupole  tom  has  its  greatest  influence  when 
die  distance  between  the  object  and  magnetic  sensor  is  less  than  the  length  of  the 
object  In  this  case,  length  of  the  105 -mm  projectile  is  42cm  and  “sensor”  to  cento 
of  spheroid  distance  is  30.25cm.  The  model  data  are  generated  at  twice  the  density 
of  the  field  data  and  it  is  not  certain  that  the  influence  of  the  octupole  tom  would  be 
observed  if  the  field  data  were  collected  at  a  greater  density.  The  Hy-Hx  component 
field  curve  is  elongated  in  the  x-direction  as  predicted  by  die  model,  however  it  does 
exhibit  some  noise  in  the  data.  In  Figure  5b  it  can  be  seen  that  the  Hz-Hx  model 
curve  approximates  the  field  curve  quite  well  when  the  projectile  has  a  90°  azimuth 
orientation  Note  there  is  no  visible  presence  of  the  octupole  term  in  the  model 
curve,  as  was  seen  when  the  projectile  was  oriented  north  (Az  =  0°).  The  Hy-Hx 
pppripnnont  plot  is  quite  interesting.  The  model  curve  indicates  there  should  be  no 
variation  in  the  y-component,  resulting  in  a  horizontal  line,  yet  the  field  curve  shows 
variation  in  the  y-component  about  three  times  greater  than  that  of  the  x-component. 
Some  of  this  variation  could  be  due  to  the  asymmetry  of  the  projectile.  However, 
such  a  large  variation  suggests  that  the  105 -mm  projectile,  which  has  not  been  fired, 
possesses  remanent  magnetization.  A  simple  test  to  check  for  remanent 
magnetisation  was  performed.  Two  north-south  profiles  were  acquired  over  the 
projectile,  one  with  the  105 -mm  oriented  pointing  north  and  the  other  south  (Figure 
6).  The  two  profiles  should  be  similar  if  no  remanent  magnetization  is  present 
These  profiles  are  radically  different,  indicating  a  strong  remanent  magnetization 
component 

An  example  of  no  remanent  magnetization  present  is  shown  in  Figure  7  for  a 
non-ordnance  steel  projectile.  The  projectile  has  a  length  of  36cm,  diameter  of  5cm, 
and  was  filled  with  sand  prior  to  firing  into  concrete.  Four  profiles  with  the 
projectile  pointing  north,  south,  east  and  west  are  plotted.  The  profile  acquired 
when  the  projectile  is  pointing  south  should  be  similar  to  the  north-pointing  profile, 
and  the  east  and  west  profiles  should  be  similar.  As  can  be  seen,  the  two  north- 
south  profiles  are  similar  in  shape  and  likewise  the  two  east-west  profiles  are 
similar.  The  impact  of  the  projectile  with  the  concrete  has  realigned  the  magnetic 
dipoles,  removing  the  presence  of  any  remanent  magnetization.  To  determine  the 
remanent  magnetization  vector,  it  is  not  sufficient  to  measure  a  single  profile  over 
the  projectile.  Barrow  et  al.  (2000)  suggest  acquiring  measurements  over  a  grid, 
with  the  projectile  located  at  the  cento  and  at  different  orientations,  and  fitting  a 
dipole  model  to  the  data  Determining  the  strength  and  direction  of  the  remanent 
magnetization  is  not  within  the  scope  of  this  report 

It  is  also  instructional  to  look  at  image  plots  of  the  individual  magnetic 
components.  Figure  8  shows  plots  of  both  model  and  field  data  for  the  total 
magnetic  field  Ht  and  three  components  (vertical  Hz,  horizontal  Hx,  Hy).  The 
model  plots  were  generated  using  the  known  field  parameters  and  105-mm 
dimensions  In  Figures  8a  and  8b,  the  upper  plots  are  for  an  azimuth  of  0°,  and  the 
lower  plots  for  an  arimnth  of  90°  (dip  0°  in  both  cases).  For  a  given  component,  the 
plots  for  Az  =  0°  and  Az  =  90°  use  the  same  shading  scale,  although  the  scales 
between  components  and  the  model  and  field  plots  differ.  This  gray  scale  uses 
black  to  indicate  magnetic  lows  and  white  for  magnetic  highs.  The  model  plots 
(Figure  8a)  have  a  similar  structure  for  both  azimuthal  orientations  of  the  spheroid. 
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although  the  magnitude  of  the  data  is  less  for  an  azimuth  of  90°.  The  similarity  is 
attributed  to  the  symmetry  of  the  spheroid,  having  an  approximately  north  pointing 
geomagnetic  field  (actual  declination  1°),  and  lag  between  the  magnetization  vector 
and  file  projectile  azimuth.  The  same  sensor-projectile  separation  was  used  to 
generate  these  data  as  was  used  to  acquire  the  field  data.  The  small  separation 
results  in  a  dual  (or  double)  dipole  observed  in  the  Hx  (north  component)  plot  for  a 
0°  azimuth.  The  two  dipoles  are  attributed  to  the  proximity  of  the  sensor  to  file 
projectile  and  its  ability  to  differentiate  both  ends  of  the  projectile.  At  a  greater 
separation  the  contours  would  appear  as  a  single  dipole.  When  the  projectile  is 
oriented  at  90°,  i.e.  long  axis  of  the  projectile  is  perpendicular  to  the  survey 
direction,  it  is  perceived  as  a  single  dipolar  signature.  In  the  Hy  (east)  component 
plots  a  dipolar  anomaly  is  sensed  on  both  the  west  and  east  sides  of  the  projectile, 
separated  by  a  zero  magnitude  contour  along  the  centerline.  When  the  sensor  is  on 
the  west  side  of  the  projectile  the  magnetic  influence  of  the  projectile  is  in  the 
positive  y-direction,  thus  giving  a  positive-negative  dipole,  whereas  when  the  sensor 
is  east  of  the  projectile  the  magnetic  influence  is  in  the  negative  y-direction, 
resulting  in  a  negative-positive  dipole. 

The  component  images  for  the  field  data  are  plotted  in  Figure  8b.  For  an 
azimuth  of  0°,  the  field  plots  are  similar  to  those  predicted  by  the  model,  although 
the  range  of  magnitude  is  greater  in  the  field  data.  Without  arty  knowledge  on  fire 
subsurface  placement  of  the  projectile  and  assuming  no  remanent  magnetization, 
some  general  statements  can  be  inferred  from  the  data.  The  classic  dipolar  response 
in  file  total  magnetic  field  plot  suggests  the  projectile  is  oriented  at  roughly  0°  or 
90°  azimuth.  The  dual  dipole  observed  in  the  Hx  component  indicates  the  projectile 
length  is  at  least  as  great  as  the  sensor-projectile  distance,  concluding  the  projectile 
is  relatively  shallow  and  of  moderate  length  or  it  is  a  larger  projectile  at  greater 
depth.  Inspection  of  north-south  and  east-west  profiles  over  the  object  may  aid  in 
determining  whether  the  projectile  is  oriented  at  an  azimuth  of  0°  or  90°.  Model 
profiles  extracted  from  the  centerlines  of  the  total  magnetic  field  (Ht)  grids  in  Figure 
8  are  plotted  in  Figure  9  along  with  the  field  curves.  Figure  9a  compares  the  model 
profiles  for  an  azimuth  of  0°  to  the  field  profiles,  and  those  in  Figure  9b  are  for  a 
90°  azimuth  As  expected,  the  model  curves  for  an  azimuth  of  0°  provide  the  better 
fit  An  attempt  to  improve  the  fit  for  the  model  Az  =  90°  by  decreasing  the 
projectile  depth  would  increase  the  magnitude  of  the  total  magnetic  field,  resulting 
in  a  better  fit  to  the  north-south  profile  and  worse  fit  for  the  east-west  profile.  Thus, 
the  general  conclusion  is  the  projectile  is  oriented  with  its  long  axis  parallel  to  the 
geomagnetic  field,  Az  =  0°.  The  above  discussion  is  based  on  the  assumption  that 
the  projectile  has  no  permanent  magnetization.  In  this  case,  the  105-mm  under 
study  has  not  been  fired  and  is  known  to  have  a  significant  remanent  magnetization 
component.  However,  a  detailed  study  by  Barrow  and  Nelson  (2000)  involving  the 
collection  and  analysis  of  spatially  dense  magnetometer  data  and  ground  truth  of 
many  of  the  anomalies,  revealed  that  ordnance  intact  after  impact  showed  little  or  no 
remanent  magnetization,  whereas  ordnance  debris  and  non-ordnance  items  were 
relatively  magnetic.  Therefore,  the  assumption  that  the  measured  magnetic  field 
over  buried  UXO  is  primarily  in  response  to  the  induced  magnetization  is  valid. 

The  data  acquired  over  the  105-mm  projectile  when  it  is  oriented  at  an  azimuth 
of  90°  (Figure  8b)  is  quite  different  than  the  model  (Figure  8a),  strongly  suggesting 
the  presence  of  remanent  magnetization.  The  aspect  ratio,  length  to  diameter,  of 
many  UXO  is  in  the  range  3  to  6.  Altshuler  (1996)  shows  that  in  this  range  the 
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tnavimiim  angle  achieved  between  the  long  axis  of  the  projectile  and  the 
geomagnetic  field  is  57  to  73  degrees,  respectively,  beyond  which  the  magnetization 
vector  moves  back  towards  the  direction  of  the  geomagnetic  field.  The  horizontally 
oriented  dipole  observed  in  the  total  magnetic  field  (Ht)  plot  clearly  indicates  a 
strong  remanent  magnetization  component.  A  prolate  spheroid  will  become 
magnetized  more  easily  along  its  long  axis  than  its  semiminor  axis  (Altshuler  1996). 
Thug  it  can  be  assumed  that  any  remanent  magnetization  present  is  likely  in  the 
direction  of  the  semimajor  axis.  Based  on  the  total  magnetic  field  image  plot,  it  can 
be  conclude  that  remanent  magnetization  is  present  and  that  the  object  is  buried  in 
an  east-west  direction.  Note  that  the  behavior  of  the  x-  and  y-component  is  opposite 
that  for  an  azimuth  of  0°,  another  indication  that  the  orientation  is  east-west.  Recall 
in  Figures  5a  and  9a  that  the  range  in  magnitude  of  the  field  data  is  greater  than  that 
predicted  by  the  model.  A  remanent  magnetization  component  aligned  with  the  long 
axis  of  the  projectile  reinforces  the  induced  magnetic  field  when  the  projectile  is 
directed  north.  Without  knowing  the  actual  orientation  of  the  105-mm  projectile, 

one  can  state  that  the  object  causing  the  anomaly  does  possess  remanent 

magnetization  but  it  is  questionable  whether  the  degree  of  remanent  magnetization 
can  be  ascertained  to  allow  determination  of  the  true  projectile  position.  If  no 
remanent  magnetization  is  assumed,  the  north-south  model  profile  for  an  azimuth  of 
90°  provides  a  relatively  good  fit  to  the  field  data,  however  there  is  no  similarity 
between  the  model  and  field  east-west  profiles  (Figure  10).  Based  on  the  findings  of 

Barrow  and  Nelson  (2000)  that  fired  UXO  exhibit  little  or  no  remanent 

magnetization,  it  may  be  concluded  that  the  anomaly  caused  by  this  projectile  is  not 
a  UXO  and  poses  no  threat  Although  this  conclusion  would  be  correct  since  this 
particular  105-mm  projectile  is  inert,  the  discovery  of  an  intact  projectile  in  the 
subsurface  on  formerly  owned  military  properly  tends  to  alarm  the  individual  that 
wrerthftH  the  projectile,  local  governments,  and  the  general  public.  Thus,  there  is 
an  urgent  and  justifiable  need  to  delineate,  differentiate,  and  identify  anomaly 
sources. 

Sensor  Levelness 

A  primary  concern  of  using  three  component  magnetic  sensors  is  the  feasibility 
of  their  use  in  a  typical  field  environment  Accurate  component  measurements  are 
dependent  on  proper  orientation  and  levelness  of  the  sensor.  To  check  the  effect  of 
improper  leveling  of  the  sensor,  controlled  magnetic  field  measurements  were 
collected  along  a  5m  north-south  profile  line  with  no  ordnance  present  using  the 
measurement  stand  pictured  in  Figure  5.  The  stand  has  less  than  a  10  incline  from 
south  to  north  over  the  5m  distance.  Even  when  the  sensor  was  leveled  using  a 
bubble  level  mounted  on  the  sensor  holder,  the  sensor  was  not  completely  level.  A 
comparison  of  average,  minimum,  and  maximum  angles  calculated  from 
accelerometer  measurements  along  the  survey  grid  during  level  and  non-level  states 
is  given  in  Table  1.  The  accelerometer  measurements  acquired  when  the  sensor  was 
level  are  only  from  one  line  on  the  grid  because  there  was  little  variation  in  the 
measurements.  When  the  sensor  was  assumed  to  be  level,  the  x-  and  y-axis  were 
still  about  1.5°  off  level.  Removing  the  average  tilt  value  when  the  sensor  is  level 
from  the  average  non-level  angle  results  in  a  compensated  average  tilt  of  3.5°  for  the 
x-axis  and  0.3°  for  the  y-axis. 
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Table  1.  Tilt  Angles  Measured  for  Level  and  Non-Level  Sensor 

Sensor  Level  (°  tilt) 

Sensor  Not  Level  (°  tilt) 

Avg 

Min 

Max 

Avg 

Min 

Max 

AVg  -  AVQLevel 

x-direction 

1.61 

1.53 

1.68 

5.12 

1.16 

7.10 

3.51 

y-direction 

1.52 

1.48 

1.56 

1.18 

-0.17 

2.58 

0.34 

The  three  magnetic  field  components  were  measured  with  the  sensor  holder 
both  level  and  not  level.  These  data  and  the  calculated  total  magnetic  field  are 
presented  in  Figure  1 1.  Variation  in  the  individual  components  when  the  sensor  is 
not  level  ranges  from  approximately  3  to  18  times  that  with  it  level.  Table  2  gives 
the  variation  observed  for  each  component  under  both  level  and  non-level 
conditions.  The  large  variation  in  the  Hy  component  level  data  is  likely  attributed 
to  the  measurement  stand  having  less  stability  in  that  direction.  Calculation  of  the 
total  magnetic  field  from  the  components  tends  to  minimize  the  non-level  effects, 
with  the  variation  in  the  non-level  data  about  three  times  greater  than  the  level  data. 
Further  study  is  needed  to  determine  if  there  is  a  relationship  that  can  be  exploited 
to  mathematically  compensate  for  a  non-level  sensor. 


Table  2.  Variation  in  Measured  Magnetic  Field  Components  for  a  Level  and  Non-Level 

Sensor 

Variation  (nT) 

Hx 

Hy 

Hz 

Ht 

Not  Level 

2580 

2080 

1215 

46 

Level 

140 

718 

77 

16 

Ratio 

18.4 

2.9 

15.8 

2.9 
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4  Summary 


This  initial  study  suggests  that  three-component  magnetic  data  may  provide 
valuable  information  that  is  not  obtainable  using  conventional  total  field 
magnetometers.  The  additional  information  will  aid  in  the  discrimination  and 
identification  ofUXO.  As  seat  in  the  examples  above,  the  vertical  z-component 
tends  to  mimic  the  total  field  measurement,  although  the  magnitude  is  less.  The 
horizontal  x-  and  y-components  provide  information  that  is  not  discemable  in  the 
total  field  data,  such  as  relative  depth  and  size  of  the  buried  object  They  also 
contribute  in  determining  the  arimuthal  orientation  and  direction  of  remanent 
magnetization. 

There  are  obstacles  to  overcome  if  three-component  magnetic  sensors  are  to  be 
used  in  practical  field  situations.  The  sensors  are  highly  sensitive  to  orientation  and 
levelness.  An  improperly  oriented  and  non-level  sensor  can  cause  a  magnitude 
increase  in  the  variation  of  the  magnetic  field  measured.  One  possible  solution  is  to 
design  an  instrument  that  would  allow  the  sensor  to  “float”  in  order  to  self-level. 
Another  solution  is  to  study  the  feasibility  of  developing  mathematical  expressions 
to  correct  for  deviations  in  sensor  orientation  and  levelness. 

The  presence  of  remanent  magnetization  in  an  object  affects  all  types  of 
magnetometers:  total  field,  gradient,  and  multi-component.  If  the  findings  of 
Barrows  and  Nelson  (2000)  can  be  applied  to  all  UXO  contaminated  sites,  then  any 
anomaly  detected  that  exhibits  obvious  remanent  magnetization  can  be  eliminated  as 
a  possible  UXO.  Further  study  is  needed  to  determine  if  in  general,  UXO  lose  their 
magnetic  history  on  impact,  and  explosive  waste  acquires  a  shock  magnetization 
when  the  projectile  explodes.  Anomalies  that  have  a  weaker  remanent 
magnetization  component,  where  it  is  not  obviously  present  from  looking  at  a  plot 
of  the  field  data,  will  adversely  affect  the  ability  to  correctly  determine  the  depth, 
azimuth  and  dip  of  an  anomaly  source.  It  is  questionable  whether  the  magnitude  of 
the  remanent  magnetization  can  be  determined  and  removed  from  the  data,  leaving 

i 

is 

remanent  magnetization  and  its  influence  on  UXO,  ordnance  waste,  man-made 
objects,  and  geologic  material  are  needed. 


only  the  indi^  magnetic  portion.  Three-component  magnetometers  do  provide  i 
greater  chance  of  identifying  and  isolating  the  remanent  magnetization  because  it 
nossible  to  see  how  the  individual  components  are  affected.  More  studies  on 
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105mm  -  Model  Data 

Incl  =  65,  Decl  =  0,  Fo  =  50000nT,  z  =  50cm 
Total  Magnetic  Field 


rotal  magnetic  field  model  data  generated  for  a  105-mm  spheroid  (length  42cm,  diameter  10.5cm).  Parameters  are:  earth's 
field  50,000nT,  inclination  65,  declination  0,  distance  between  sensor  and  center  of  projectile  50cm,  magnetic 
permeability  1000.  (a)  spheroid  azimuth  varies,  spheroid  dip  =  0.  (b)  spheroid  dip  varies,  spheroid  azimuth  =  0. 


(b)  Azimuth  = 


105-mm  Model  Data  w /  Spherical  Model  Fit 
105-mm:  Dia  =  10.5cm,  L  =  42cm;  Sphere:  Dia  =  15cm 


(iu)  Play  0ijau6ew  |e»oj. 
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Figure  2.  Spherical  model  fit  to  105-mm  spheroid  model  profiles  in  Figure  1.  Sphere  diameter  15cm.  (a)  spheroid 
azimuth  varies,  spheroid  dip  -  0.  (b)  spheroid  dip  varies,  spheroid  azimuth  =  0. 


105-mm  Model  Data 

Incl  =  65,  Decl  =  0,  Fo  =  50000nT,  z=  50  cm 


(b)  spheroid  dip  varies,  apheroid  azimuth  = 


(bj  magnetometer  sensor  in  holder. 


Figure  4.  Concluded 


105-mm  Field  Data  w /  Spheroid  Model  Fit 
Incl  =  62,  Decl  =  0,  Dip  =  0,  z  =  30.26cm 


Figure  5.  Magnetic  component  plots  for  105-mm  field  data  and  spheroid  model.  Parameters  are:  earth's  field  50,200nT, 

inclination  62,  declination  0,  distance  between  sensor  and  center  of  projectile  30.25cm,  magnetic  permeability  1000. 
(a)  projectile  azimuth  =  0,  dip  =  0.  (b)  projectile  azimuth  =  90,  dip  =  0. 


105-mm  Projectile 


Figure  6.  Profiles  acquired  over  105-mm  projectile  with  projectile  pointing  north  and 
south  demonstrating  presence  of  remanent  magnetization 


Steel  Projectile,  Fired 


Figure  7.  Example  of  no  remanent  magnetization  as  demonstrated  by  profiles  acquired 
over  a  fired  steel  projectile 


Image  plots  of  the  total  magnetic  field  and  individual  components  Hx,  Hy,  Hz  for  (a)  105-mm  spheroid 
model  and  (b)  105-mm  projectile  field  data 


1 05-mm  Projectile  Field  Data 


105-mm  projectile,  Az  =  0 
Spheroid  model,  Az  =  0 

Profiles  Along  Centerline  of  Ht  Plots  in  Figure  8 
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Figure  9.  Comparison  of  total  magnetic  field  profiles  for  model  and  field  data  in  Figure  8.  (a)  105-mm  projectile  azimuth  =  0, 
spheroid  azimuth  -  0.  (b)  105-mm  projectile  azimuth  =  0,  spheroid  azimuth  =  90. 
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Figure  10.  Comparison  of  total  magnetic  field  profiles  of  data  in  Figure  8  for  both  model  spheroid  and  105-mm  projectile 
having  an  azimuth  of  90 
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